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ABSTRACT 

The  stability  has  been  investigated  of  the  unsteady  flow  past  an  infinite  flat  plate 
when  it  is  moved  impulsively  from  rest,  in  its  own  plane.  For  small  times  the  in¬ 
stantaneous  stability  of  the  flow  depends  on  the  linearised  equations  of  motion  which 
reduce  in  this  problem  to  the  Orr-Sommerfeld  equation.  It  is  known  that  the  flow  for 
certain  values  of  Reynolds  number,  frequency  and  wavenumber  is  unstable  to  Tollmien- 
Schlichting  waves,  as  in  the  case  of  the  Blasius  boundary  layer  flow  past  a  flat  plate. 
With  increase  in  time,  the  unstable  waves  only  undergo  growth  for  a  finite  time  inter¬ 
val,  and  this  growth  rate  is  itself  a  function  of  time.  The  influence  of  finite  amplitude 
effects  is  studied  by  solving  the  full  Navier-Stokes  equations.  It  is  found  that  the 
stability  characteristics  are  markedly  changed  both  by  the  consideration  of  the  time 
evolution  of  the  flow,  and  by  the  introduction  of  finite  amplitude  effects. 

*  Research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  contract 
No.  NASl-19480  while  the  authors  were  in  residence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681 
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§1  Introduction 


In  the  flow  over  a  flat  plate  at  zero  pressure  gradient  the  stability  of  the  boundary 
layer  is  not  only  dependent  on  its  instantaneous  thickness,  but  also  on  its  rate  of 
growth.  This  problem,  based  on  the  Blasius  mean  flow  has  been  studied  extensively. 
However,  experiment  shows  that  the  stability  characteristics  based  on  the  local  growth 
of  the  boundary  layer  provide  a  good  approximate  model.  This  weis  confirmed  in  the 
triple-deck  structure  calculations  of  Smith  (1979a)  on  the  influence  of  the  rate  of 
growth  of  near  neutral  unstable  Tollmien-Schlichting  (TS)  modes  at  high  Reynolds 
numbers.  The  weakly  nonlinear  evolution  of  TS  modes  within  a  Blasius  layer  was 
studied  by  Smith  (1979b)  and  Hall  &  Smith  (1984).  In  the  latter  of  these  studies,  it 
was  demonstrated  how  a  linear  mode  could  evolve  smoothly  into  a  finite  amplitude 
perturbation.  In  Otto  (1994),  this  was  shown  to  be  feasible  within  the  unsteady 
problem.  Smith  also  found  that  it  was  not  possible  to  determine  the  effects  of  boundary 
layer  spatial  growth  on  the  stability  of  the  modes  since  different  conclusions  were 
foimd  depending  on  which  flow  quantity  was  used  to  express  the  finite  amplitude  of 
a  mode.  Caster  (1974)  and  Eagles  &  Weissman  (1975)  have  edso  contributed  to  our 
imderstzmding  of  the  Blasius  problem.  In  this  flow,  Tollmien-Schlichting  waves  grow 
.spatially  in  the  flow  direction  and  the  resulting  equations  in  this  coordinate  are  elliptic, 
and  are  computationally  expensive  to  solve.  The  approach  of  Bertolotti  (1991),  using 
the  Parabolised  Stability  Equations  (PSE),  in  which  the  elliptic  equations  are  made 
parabolic  on  the  assumption  of  slowly  evolving  wave  characteristics  and  the  changes  in 
structure  may  be  locally  zissumed  linear,  has  been  one  of  the  more  successful  methods 
to  reduce  the  computational  cost.  This  approach  allows  feir  larger  steps  to  be  taken  in 
the  evolutionary  coordinate. 

In  the  present  work,  the  related  temporal  problem  is  considered.  Otto  (1994)  found 
rea.sonable  agreement  between  triple-deck  and  Orr-Sommerfeld  results  in  regimes  of 
common  validity,  for  the  unsteady  flow  problem  of  the  flat  plate  moved  impulsively 
from  rest.  Here,  we  discuss  both  the  linear  and  nonlinear  problems,  and  since  our 
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problem  is  already  parabolic  in  the  evolutionary  coordinate,  we  use  the  PSE  method¬ 
ology  to  study  the  effects  of  finite  amplitude  and  to  reduce  the  number  of  time  steps 
per  period  required  for  accuracy,  which  it  does  very  effectively.  We  find  substajitial 
differences  between  the  results  of  linear  theory  based  on  the  Orr-Sommerfeld  equations 
and  our  results  for  finite  amplitude  waves. 

In  Section  2,  we  derive  equations  for  the  linear  and  nonlinear  evolution  of  Tollmien- 
Schlichting  waves  with  a  brief  summary  of  our  period  fitting  technique.  Section  3 
contains  details  of  the  numerical  methods  used  to  solve  the  equations.  In  Section  4, 
we  discuss  the  results  and  finally  in  Section  5,  we  draw  some  conclusions. 
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§2  Formulation 

The  physical  problem  consists  of  an  infinite  flat  plate  in  two  dimensions  which  is 
at  rest  within  a  fluid  of  density  po  and  pressure  po-  At  time  t  =  0,  the  plate  is  moved 
impulsively  from  rest  with  the  speed  Uq,  in  the  direction  parallel  to  its  alignment  (see 
Figure  2.1). 


The  governing  equations  for  this  problem  are  the  equation  of  continuity  and  the 
Navier-Stokes  equations 


du  ^  _  n 

dx  ^  dy 

du  du  du  dp  1  _2 

dt^  dx^  dy  dx^Re 

dv  dv  dv  dp  1  _2 

dt^  dx^  dy  dy^Re 


(2.1) 


where  each  variable  is  dimensionless,  with  scaling  parameters  po,  Uq,  po,  and  L  repre¬ 
senting  the  ambient  density,  plate  velocity,  viscosity  and  a  unit  of  length  in  the  flow 
direction,  respectively.  The  Reynolds  number  per  unit  length  of  plate,  Re  is  defined 
as 


Re  = 


LUopo 

fio 


(2,2) 


and  is  the  two-dimensional  Laplacian  operator 

= - 1 - 

9*2  ^  ay2 


(2.3) 


We  note  that  the  pressure  has  been  non-dimensionalised  by  PqUq  and  time  by 


VSL 
L  • 
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In  the  case  of  aai  infinite  plate  impulsively  set  in  motion,  the  mean  solution,  i.e. 
the  undisturbed  flow,  is  given  in  dimensionless  terms  by  the  reduced  set  of  equations 


du  1  d^ii 

dt  Re  dy^ 

u  =  0 


(2.4) 


P  = 


Po 

p«ui 


with  initial  and  boimdary  conditions 

u(0,()  =  |“  \ 


=  0 
>0 

u{Y,t)  — »  0  K  — oo 


(2.5) 


The  mean  streamwise  velocity  component  is  therefore  given  exactly  by 


(2.6) 


Given  the  form  of  (2.6),  it  is  appropriate  to  choose  physical  variables  X  and  Y  that 
absorb  the  \/Sc  dependence, 

X=xVl^ 


Y  =  yy/R^ 

With  this  transformation,  the  governing  set  of  equations  become 


(2.7) 


du 

1 

dv 

_ n 

dX 

dY 

=  u 

1 

du 

du 

du 

dp 

1 

dt 

dY 

dX 

y/Re 

1 

dv 

dv 

dv 

dp 

1  , 

y/Re 

dt 

dY 

dY 

4- 

y/R/ 

(2.8) 


and  the  mean  flow  is  now  given  by 


«(r,<)  =  l-erf(^). 


(2.9) 


The  Izwjk  of  a  physical  boundary  condition  for  pressure  provides  a  computational 
diflBctilty  for  equations  (2.8)  cast  in  terms  of  the  primitive  variables.  We  avoid  this 
problem  by  reducing  this  set  of  equations  to  a  single  equation  for  the  normal  velocity 
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component.  In  this  section,  we  will  tliscuss  the  analysis  for  both  linear  and  nonlinear 
disturbances.  In  addition,  we  show  how  to  generate  an  appropriate  initial  condition 
to  begin  the  stability  calculation,  and  introduce  a  period-fitting  technique  intended  to 
reduce  the  total  computational  effort. 

2.1  Linear  Disturbance  Equations 

We  let  u,  V  and  p  represent  disturbances  from  the  mean  state: 


u{X,Y,t)  =  u{Y,t)  +  u{X,Y,t), 
v{X,Y,t)  =  v{X,Yj), 
p(X,Y,f\  =  -^+p{X,Yj). 


(2.10) 


After  substituting  equations  (2.10)  into  equations  (2.8)  and  linearizing  the  resulting 
expressions,  the  perturbation  pressure,  p,  and  the  strezunwise  velocity  component,  u, 
may  be  eliminated,  giving  the  following  single  equation  for  v,  the  transverse  compo¬ 
nent: 

f/  1  d  _  _  d  1  ,  ,92,  d^u  d]  ^  ^ 

\y/R^dt^''dx  9^-2  +  +  9^2 )  dY'^dxy  (  •  ) 

The  physical  problem  is  spatially  homogeneous  in  the  X-direction,  zind  it  is  legitimate 
to  assume  a  periodic  X-variation,  giving 


v(X,F,t)  =  u(F,t)e‘ 


(2.12) 


which  reduces  (2.11)  to 


^y/Wtdt  "“^^9y2  ""^'^"9y2 202^^2+®")  ’'  0.  (2.13) 


When  a  periodic  time  dependence,  e~*‘*'‘,  is  assumed,  and  noting  the  mean  velocity 
is  a  function  of  Y  and  t  only,  equation  (2.13)  reduces  to  the  familiar  Orr-Sommerfeld 
equation  with  similar  results.  In  our  case,  we  consider  the  solution  to  equation  (2.13) 
at  zmy  fixed  time  and  then  the  results  show  the  local  stability  characteristics  of  the 


flow  for  the  linear  problem.  For  this  time- dependent  problem,  the  boundary  conditions 
are 


and  since  ^  =  0  at  F  —  0, 


v(0,t)  =  0. 


dv 

W 


(0,i)  =  0. 


(2.14) 


(2.15) 


In  the  farfield,  as  F  oo,  the  solution  depends  on  e  giving  for  large  F, 


+  av(YJ)  =  0 

5^-a^«(y,()  =  0 


OO  . 


(2.16) 


2.2  Nonlinear  Disturbance  Equations 

Substituting  (2.10)  into  (2.8),  retaining  the  nonlinear  terms,  and  eliminating  the 
perturbation  pressure,  we  find  the  following  for  v: 

(Id  _  d  ^  d'^  ^  d'^u  d 

^V^dt'^'^dX  dY^dXy~ 

( ^  d  ~  d  \  /  d'^  d'^  \  ,  dv  /  \  -  dv  f  d'^  d'^  \_ 

~  v“lx  ■'■''aF/  vajfs  ajn) "  W^r  ax  vax*  W^r' 

(2.17) 

The  could  be  eliminated  using  the  continuity  equation.  However,  since  the 
term  cannot,  we  leave  both  terms  in  for  ease  of  presentation. 

We  now  express  the  disturbance  quantities  as  a  finite  sum  of  modes  with  respect 
to  the  wavenumber  alpha: 

M 

i{X,Y,t)=  Y, 

j=~M 

M 

u{X,Y,t)=  Y,  «;(!'.  (2.18) 

j=~M 

M 

p(x,y()=  Y  P;(n 

j=-M 

The  continuity  equation  allows  for  elimination  of  all  the  u/s,  except  uq,  in  favour  of 
the  tij’s  and  requires  that  VQ{Y,t)  =  0.  The  resulting  equation  after  making  these 
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substitutions  in  (2.17)  and  neglecting  all  products  which  combine  to  give  a  mode 
beyond  M  is  given  by 


/  i  d 

\\/Re 


1  ,  •>  9\\  /  9  9\  d^u] 


where 


E.,  k  -  ifVj-k\  ,  J 

Jes~-  M  V.  *f  >f 


k  ,  dvj-k 


j  d^Vk  dvj-k  j  d^vk  . 


k  dY^  dY  k  dY^ 


■  -2  2-- 
+  -  J  a  uoVj 


d'^uo  . 


(2.20) 


for  j  =  — M,  —M  +  1,  —M  +  2, ...,  —1,  and 


Ni)  = 


2^  ■'  \dY  ’  i-k  dY  j-kdY^  dY  k  dY> 


k  =  -M+j 
kjtO,} 


,  .  /.  d^Vj  -2  2..  d^Uo  . 


(2.21) 


for  J  =  1,2, 3,...,  M. 

The  boundary  conditions  for  the  Vj's  axe  analagous  to  those  for  the  linear  case: 


0,(0,t)  =  0,  ||-(0,t)  =  0, 

-%y’"  +jocvjiY,t)  =  0 
-  pa^v^(Y,t)  =  0 


(2.22) 


y  — >■  oc . 


The  variable  uo{Y,t)  has  no  associated  wavy  part,  and  is  effectively  a  correction  to 
the  mean  solution  due  to  the  presence  of  the  finite  amplitude  waves.  An  equation  for 
it  may  be  found  by  looking  at  the  non-wavy  component  of  the  x-momentum  equation, 
after  substituting  in  (2.18): 


duo  d^iiQ 
~dt  ~  dY^ 


Z.  kadY^''~'‘' 


(2.23) 
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The  boundary  conditions  for  uq  consist  of  a  no-slip  condition  at  tlie  wall  and  an 
asymptotic  condition  far  away  from  the  wall: 

tio(0,<)  =  0,  =0  Y^oo.  (2.24) 

The  system  of  equations  (2.19)  through  (2.24)  define  the  nonlinear  problem. 

2.3  Initial  Condition 

In  order  to  solve  either  the  linear  or  nonlinear  problem  for  some  time  t  =  to  >  0, 
we  need  to  define  an  initial  perturbation  which  is  at  least  an  approximate  solution 
of  the  linearised  equations.  To  this  end,  we  can  freeze  the  flow  at  a  given  time  and 
determine  the  local  stability  characteristics  of  the  flow  from  an  Orr-Sommerfeld  (OS) 
analysis.  In  particular,  for  given  values  of  a  and  v/^,  we  can  determine  a  function 
u(F)  and  a  frequency  uj,  such  that  an  initial  disturbance  of  the  form  u(y)e‘"^  will 
behave  according  to  for  all  time  for  the  frozen  flow. 

However,  since  the  OS  equation  does  not  allow  for  the  evolution  of  the  undisturbed 
flow  with  time,  using  an  initial  condition  constructed  in  this  way  generates  spurious 
oscillations  in  the  solution.  We  overcome  this  problem  by  finding  an  initial  condition, 
which  is  a  slight  perturbation  of  the  OS  eigenfunction  u(F),  that  is  sensitive  to  the 
changing  mean  flow.  Following  the  approach  of  Bertolotti,  we  let 

v{Y,  t)  =  v'{Y,  t)e~'  fo  (2.25) 


where  it  is  assumed  that  the  time  dependence  of  v'  is  slowly  varying.  Substituting  this 
expression  into  the  linear  stability  equation  (2.13)  gives 


( 


y/Re 


av'(Y,t) 

at 


(2.26) 
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We  now  make  a  order  Taylor  Series  ai)proximation  about  to- 


t  S3  to  +  T  , 


6  Si  9{to)  +  —  ^0  +  , 

t;'(y,t)  «  r;'(y, to)  +  r^(y  to)  =  v'^iV)  +  ri>;(r) , 


at 


dv' 

at 


(r,io)  =  v{(r), 


(2.27) 


u(y  t)  %  M(y  to)  +  r|^(y  to), 


^  aiu 


ays  "  ~  3v^{y  ^o)  +  to) . 

We  note  that  the  approximation  for  does  not  include  a  ^j^(y  to)  term;  the 

slow  variation  of  v'  with  t  allows  us  to  ignore  the  second  derivative.  Collecting  like 
powers  of  r,  namely  r®  and  t\  we  have  the  following  equations,  respectively: 

LvUy)  =  -i(^-a^)v',(y) 

-  “")*’;(r)  -  a5gf,(r,foK(r) 

where  L  is  the  Orr-Sommerfeld  operator 


(2.28) 


^  +*Jy7(y  to)+ 2q;^  gy2  +0;'*^.  (2.29) 

The  boundary  conditions  for  Vg  and  Uj  are  given  by  (2.14),  (2.15),  and  (2.16). 

With  the  additional  constraint  that  Vg  arid  v\  are  orthogonal, 

(t;;)*v;dy  =  0,  (2.30) 

equations  (2.28)  provide  an  eigenvalue  problem  which,  when  solved,  gives  an  initial 
profile  Vg{Y)  and  frequency  9q  which,  to  first  order,  take  into  account  the  evolution  of 
the  undisturbed  flow. 

To  further  improve  the  quality  of  the  initial  condition,  a  higher  order  version  of 
the  above  can  be  performed.  In  this  case,  we  use  a  second  order  Taylor  Series  approx¬ 
imation,  and  include  the  second  derivative  of  v'  with  respect  to  t,  but  ignore  the  third 
derivative.  In  addition  to  Vg  and  uj,  the  resulting  eigenvalue  problem  contains 

unknowns  ^2  and  which  represent  the  second-order  Taylor  Series  terms.  Extending 
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the  above  analysis  to  2"*^  order  gives  a  linear  equation  for  and  an  orthogonality 
condition,  analagous  to  (2.30). 

2.4  Period  Fitting 

As  mentioned  in  the  previous  section,  we  can  .separate  v(V,t),  for  the  linear  case, 
into  a  time  periodic  part,  and  a  function  v'{Y,t)  that  varies  slowly  with  time.  In  this 
case,  we  will  write  v  as 

fi(F,t)  =  u'(F,i)e-“^'  (2.31) 


where  u)  is  taJcen  to  be  a  real  quantity;  the  growth  or  decay  is  represented  by  the 
t- variation  of  v'.  To  adequately  resolve  the  oscillatory  part  requires  many  time  steps 
per  period  of  oscillation.  However,  the  temporal  resolution  requirements  of  v'  are 
far  less  restrictive.  Consequently,  it  is  more  efficient  to  remove  the  periodic  paxt  of 
the  solution  from  the  equations  and  compute  only  the  slowly  varying  function  v'. 
Substituting  (2.31)  into  the  linear  disturbance  equation  (2.13),  we  find  the  following 
equation  for  v' 


(2.32) 


Since  the  extent  of  the  flow  increases  as  y/i,  the  value  of  u  changes  with  time. 
Therefore,  to  implement  (2.32),  we  must  be  able  to  estimate  u;  at  any  given  time.  We 
define  the  energy  E{Y,  t)  in  the  following  way: 


roo  AOO 

E{Y,t)=  /  v\Y,t)dY=  /  u'2(F,t)e-2‘‘^‘dF.  (2.33) 

Jo  Jo 

Ignoring  the  slow  changes  of  v'  and  w  with  time,  we  approximate  a;  by  —  ^  . 

However,  we  wish  now  to  work  only  with  v'.  Initially,  we  use  a  value  of  u)  in  (2.32) 
given  by  the  OS  mialysis  or  by  the  method  described  in  Section  2.3  above.  At  each 
time  step,  we  can  correct  this  value  by  recognizing  that  the  computed  value  of  v' 
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will  contain  the  slowly  changing  envelope  plus  an  oscillatory  part  with  a  very  small 
frequency  that  represents  the  change  in  u;,  i.e. 


E'{Y,t)=  /  v^\Y,t)dY, 
Jo 


^  —  ' 


1  dE' 


2iE>  at  ’ 

+  Au;. 


(2.34) 


For  the  nonlinear  case,  the  procedure  is  very  similar.  Given  that  we  have 


aj  =  ja 


(2.35) 


we  make  the  assumption 


Uj  =  jto. 


(2.36) 


That  is,  we  compute  u)\ ,  associated  with  the  fundamental  mode,  by  the  method  shown 
in  equation  (2.34).  Then  all  ujj  for  higher  order  modes  are  given  by  (2.36).  We  note 
that  any  errors  incurred  by  this  approximation  are  at  least  partially  corrected  by  the 
the  time  dependence  of  u'  . 

The  nonlinear  disturbance  equations  for  the  ujs  now  become 


-;V) 


9v'(V,t) 

at 


+j'‘a'‘)t>'{K,()  =  Nu' 


(2.37) 

where  Nv'  has  the  same  form  as  given  in  (2.20)  except  that  v  is  replaced  everywhere 
by  v'. 
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§3  Numerical  Methods 


In  this  section,  we  will  provide  some  of  the  pertinent  details  of  the  numerical 
solution  of  both  the  linear  and  nonlinear  problems.  We  begin  with  a  brief  discussion  of 
the  grid  and  spectral  discretization,  and  follow  with  descriptions  of  the  initial  condition 
generation,  linear  simulations  and  nonlinear  simulations. 

3.1  Discretization 

The  mean  solution  is  given  by  equation  (2.9)  which  is  repeated  here: 

i2(y,t)  =  l-erf(^).  (3.1) 

Although,  u  is  technically  finite  for  any  finite  Y  and  f,  it  decays  quickly  away  from  the 
plate,  and  from  a  computational  standpoint,  it  is  appropriate  to  define  the  extent  of 
the  undisturbed  flow  to  be  the  Y  value  at  which  the  undisturbed  velocity  is  1%  of  the 
plate  speed.  Since  typically  the  OS  eigensolutions  decay  slowly  relative  to  the  decay  of 
the  undisturbed  flow  solution,  it  is  necessary  to  choose  a  physical  domain  that  extends 
well  beyond  this  limit.  Defining  the  edge  of  the  undisturbed  flow  at  t  =  to  to  be 

r.d  5  =  1  -  erf  (5^)  =  O.OI,  (3.2) 

we  have  Y^d.  =  1.821(2-^/^).  We  choose  ein  upper  bound  for  Y  to  be 

Tinax  =  40v^,  (3.3) 

i.e.  our  domain  extends  to  about  10  times  the  initial  undisturbed  flow  extent. 

We  employ  a  Chebyshev  collocation  discretization  using  iV  =  80  polynomials.  Al¬ 
though  the  disturbances  do  not  require  any  local  resolution  because,  once  introduced, 
they  extend  throughout  the  domain,  the  mean  solution  decays  rapidly  away  from  the 
wadi.  Consequently,  we  use  a  mapping  which  maps  0  to  yinax  onto  [-1,1]  as  required 
by  the  spectral  discretization,  but  does  so  in  such  a  way  that  points  are  concentrated 
near  the  lower  boundary.  The  computational  coordinate  tj  is  defined  as  follows: 

1  =  W  (3-4) 
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where  b  and  c  are  chosen  so  that  F(l)  =  Vlnax  and  F(0)  =  Kmid  The  value  of  Kmid  is 
the  physical  coordinate  that  bisects  the  grid.  We  taJse  it  to  be  IS-y/to,  or  about  4  times 
the  initial  undisturbed  flow  extent.  This  allows  for  the  undisturbed  flow  to  grow  in 
extent  during  the  simulation,  and  keeps  the  stretching  from  being  too  severe.  With 
these  stipulations  we  have 


Figure  3. 1  Representation  of  Stretched  Spectral  Grid 


3.2  Initial  Conditions 

For  both  the  lineax  and  nonlinear  simulations,  the  disturbances  are  introduced  at 
dimensionless  time  to  =  0.25.  We  generate  three  sets  of  initial  conditions  at  this  fo- 
The  first  corresponds  to  the  OS  eigenfunction,  and  the  remaining  two  eire  the  and 
2”**  order  corrections,  respectively,  of  §2.3. 

The  continuous  Orr-Sommerfeld  eigenvalue  problem  is  defined  by 

Lu' (y)  =i  0  (3.6) 


and  the  boundary  conditions 

>’o(0)=0,  1^(0)  =  0, 

2^-»X(l')  =  0  ’ 

where  L  is  the  Orr-Sommerfeld  operator  given  by  equation  (2.29).  The  corresponding 
discrete  problem  replaces  ^  with  the  discrete  operator  (^),,  where  given  a  discrete 
function  <f>  defined  on  the  nodes  I  =  0, 1,2, ...,  N, 

(^)/  =  (If)/  X!  (3-^) 
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where,  once  again,  t]  G  [—1,1]  is  the  computational  coordinate.  The  are  the 

coefficients  of  the  Chebyshev  1*‘  derivative  matrix.  The  system  is  then  defined  by 
applying  the  discrete  OS  equation  at  the  interior  nodes  2, 3, 4,  ...,N—2,  and  the  discrete 
wall  eind  asymptotic  boundary  conditions  at  t;  =  —  1  (F  =  0)  and  r/  =  1  (F  =  F^ax) 
respectively. 

For  a  given  a  and  y/Re,  the  solution  of  this  eigenvalue  problem  provides  a  frequency 
LJ  and  an  initial  profile  u®,/  =  0,1,2, ...,1V,  which  may  be  used  to  start  the  linear 
simulation. 

In  the  nonlinear  case,  this  profile  may  be  used  as  the  starting  value  for  the  fu 
mental  mode.  In  these  simulations,  we  set  the  higher  order  modes  initially  to  0, 
the  —1  mode  to  the  complex  conjugate  of  the  fundamental. 

The  discrete  analogues  to  the  1*‘  and  2”**  order  correction  problems  of  §2.3  follow 
directly  from  above.  The  resulting  linear  systems  are  also  eigenvalue  problems,  the 
solutions  of  which  provide  a  frequency  and  initial  profile. 


3.3  Linear  Simulations 

The  continuous  linear  problem  is  defined  by  equation  (2.13)  with  boundary  condi¬ 
tions  (2.14),  (2.15),  and  (2.16)  and  initial  condition  provided  as  in  the  the  above  sec¬ 
tion.  We  employ  the  discrete  form  of  these  equations  along  with  the  Crank-Nicolson 
time  differencing,  given  by 


h  - 


_ + 

At  ~  2  V  at  ^  at  )' 
which  results  in  the  following  linear  set  of  equations  for 


[( 


At 


2«\/S?' 


and 


(3.9) 


)((*)'-« 

0  - 

(3.10) 

-ha")|vr  /  =  2,3,4,...,iV-2, 

=  0, 

J 

(^).  =  0. 

9;+’  =  o,  ( 

=0. 
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The  value  of  N  for  zdl  simulations  is  chosen  to  be  ^  =  80. 


To  perform  the  simulation  with  period  fitting,  we  implement  the  discrete  form 
of  (2.32)  along  with  the  Crank- Nicolson  time  differencing  (3.9)  and  the  appropriate 
boundary  conditions.  This  linear  set  of  equations  for  (v')"*^^ ,  /  =  0, 1, 2, ...,  N  is  given 

by 


=  (3.11) 

-  iiTfe  +  «‘)]  (■’')”  '  =  2. 3. 4. ...,  W  -  2 , 

(v'To*' = 0,  ^(v'rr' = 0, 

+  -)(«%*'  =  0.  ((*)’  -  =  0. 


3.4  Nonlinear  Simulations 

The  continuous  equations  for  the  nonlinear  case  are  given  by  equations  (2.19) 
through  (2.24).  The  time  differencing  will  once  again  be  of  the  form  (3.9).  The  result¬ 
ing  discrete  set  of  equations  is  nonlinem:,  requiring  iterative  solution.  Consequently, 
we  will  verify  results  by  also  performing  a  simulation  which  uses  Crank-Nicolson  for 
the  linear  left-hand  side  of  equation  (2.19)  and  explicit  2”**  order  Adams-Bashforth  for 
the  nonlinear  terms  on  the  right-hand  side. 

For  the  explicit  treatment  of  the  nonlinear  terms,  we  have  the  following  discrete 
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analogue  to  equations  (2.19)  through  (2.21): 


(7k  - 

+  iafe((i^)'-2.v(^)Vjv)]i”r 

-5:efe((w)*-W(^)“  +  ;V)]i’:, 

+  ¥(3nc';,-n«;7‘) 


for  /  =  2, 3, ...,  TV  —  2,  where 


r  o/^^ki  k  Di}".,N  j  D^u^.DC^.. 

Mr.H  ■l--2<  k,l  ~n  J-fc,/  \  ,  J _ kj  J~k,l 

*2^  (  jjy  j  /■’'j-jfcDy2  Dy 


-7  0  UoV>,/  - 


«/<)’ 


(3.13) 


for  j  =  —M,  —M  +  1,  —M  +  2, ...,  —1,  and 


N.,,,  = ,  ^  Jta  -  —,,,-^y+  —  p-y-r  -py- 

L 


*#0,j 


7  z 


fe  DF^ 


D^v? 


-•n  .  ■  /-  -2  2-  -n  ^ 

-7  a  uovjj  - 


(3.14) 

for  7  =  1,2, 3,...,  M.  In  these  expressions,  the  superscript  represents  the  time  level, 
the  first  subscript  indicates  the  mode,  and  the  second  subscript  refers  to  the  grid. 
Similarly,  the  discrete  values  of  the  mean  flow  correction  are  governed  by 


(1  -  =  (1  +  x{w)')“o", 

M  ,  ,  , 

_  ^  f  3 -n  P  ”*,1  _  l.-in— 1  ^  ^*,1  ^ 

2a  2-/  Dyl  fc^(-*,Z)  bV*  )  ' 


(3.15) 


The  boundary  conditions  for  each  mode  j  €  [— M,  —  M  +  axe  given  by  the 

discrete  forms  of  (2.22)  for 


(^)« + -KV  =  0. 

and  (2.24)  for  the  mean  flow  correction: 


n  +  I 


^1,0 


=  0, 


(^)o  =  0^ 


=  0, 


(3.16) 


^n+l  _  Q  /  D“g'*'*  'I  Q 

“0,0  —  \  DY  f  N  ^ 


(3.17) 


The  above  equations  constitute  a  linear  set  of  equations  for  the  unknowns  j  = 

— M,  —M  +  1, ...,  M,  /  =  0, 1, 2, N,  where,  again,  N  is  taken  to  be  80.  We  use  this 
simulation  only  as  a  check.  The  period-fitting  technique  discussed  in  §2.4  allows  for 
a  drastic  increase  in  time  step.  However,  stability  considerations  will  not  allow  full 
advantage  of  this  increase  to  be  taken  when  the  nonlinear  terms  axe  treated  explicitly. 
The  majority  of  the  resiilts  presented  in  the  following  Results  section  were  obtained 
using  period  fitting  and  Crank-Nicolson  differencing  for  the  nonlinear  terms,  resulting 
in  a  fully  nonlinear  set  of  algebraic  equations  to  be  solved  at  each  time  step.  In  this 
case,  (3.12)  is  replaced  by  the  discrete  finalogue  of  (2.37) 


[(^fe  +  ~  >““"*'))  (( w)*  -  A' 

(3.18) 

i-')h 

+  f(N(»')r+N(v')r), 


and  (3.15)  is  replaced  by 


M 

i\/ReAt  I 

2^;  2^  k\y{-k,l)  DY^  DV^  )  • 


(3.19) 


k  =  -M 
***0 
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To  solve  the  nonlinear  set,  at  each  time  step  we  begin  by  making  a  guess  for  the 
solution  which  is  a  cubic  extrapolation  using  the  previous  4  time  steps.  Using  this 
guess,  we  compute  the  nonlineeir  terms  on  the  right-hand  sides  of  equations  (3.18)  and 
(3.19),  and  then  solve  the  resulting  linear  systems  for  new  values  of  which  serve 
as  the  new  guesses  in  the  iteration.  The  process  is  repeated  until  convergence. 

To  utilise  period-fitting,  we  employ  the  discrete  analogue  of  equation  (2.37).  The 
algorithm  for  solving  the  resulting  nonlinear  algebraic  set  of  equations  at  each  time 
step  is  then  identical  to  that  described  above. 
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§4  Results 


In  tliis  section,  we  present  results  for  both  the  linear  and  nonlinear  simulations. 
We  begin,  however,  by  presenting  the  initial  conditions  obtained  by  Orr-Sommerfeld 
analysis,  and  by  the  and  2’"^  order  corrections. 

4.1  Initial  Conditions 

The  discrete  version  of  the  Orr-Sommerfeld  equation  (3.6)  and  corresponding  bound- 
zury  conditions  (3.7)  is  solved  by  fixing  values  of  a  and  y/Re  and  solving  the  eigenvalue 
problem  for  the  frequency  u)  and  eigenfimction  Vq.  The  value  of  u;  is  generally  complex 
with  a  positive  imaginary  part  representing  growth  with  time  and  a  negative  imagi¬ 
nary  part  indicating  temporal  decay.  For  particular  pairs  of  a  and  \/^  a  purely  real 
u)  results.  Figure  (4.1)  shows  the  neutral  Orr-Sommerfeld  curve  for  time  to  =  0.25. 
The  results  are  identical  to  those  given  in  Otto  (1994).  We  present  this  figure  for  later 
reference. 


Reg 

Figure  4.1  Neutrol  Orr-Sommerfeld  curve  showing  the  pairs 
of  Reynolds  number  and  streamwise  wave  num¬ 
ber  which  give  purely  real  u 
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The  reference  values  of  a  and  \/Re  we  have  chosen  for  our  simulations  are 


a  —  0.4 
=  3000 


(4.1) 


which  defines  a  point  just  to  the  left  of  the  neutral  curve  in  figure  (4.1).  For  these 
values,  the  corresponding  uj  resulting  from  the  OS  solution  and  the  1*‘  and  2"“*  order 
corrections  axe  given  by 

Orr-Sommerfeld:  0.27589729  -  0.00004782i 

1*‘  Order  Correction  0.27585418  +  0.00028591i 

2”''  Order  Correction  0.27584836  +  0.00028872z 

and  the  resulting  eigenfunctions,  which  may  be  used  as  initial  conditions  to  the  linear 
and  nonlinear  simulations,  are  presented  in  figures  (4.2a)  and  (4.2b).  Figure  (4.2a) 
shows  the  real  part  of  the  initial  profile,  as  generated  by  the  OS  analysis,  which  we 
call  In  figure  (4.2b),  two  curves  are  given.  The  dashed  curve  shows  the  deviation 
of  the  real  part  of  the  initial  profile  generated  by  the  1**  order  correction,  Vq,  from 
that  generated  by  the  OS  analysis,  Vq^.  The  dotted  curve  in  this  figure  shows  the 
deviation  of  the  real  part  of  the  initial  profile  generated  by  the  2"'*  order  correction, 
Vq,  from  that  given  by  the  1'*  order  correction.  To  plot  both  of  these  on  the  same 
graph,  we  scale  the  two  curves  appropriately,  so  the  actual  quzintities  seen  on  this 
figure  are  |10® |  and  llO"*-? — “-j  vs.  Y,  where  the  denominators  represent  the 
peaks  of  these  curves.  This  figure  shows  that  the  1**  order  correction  represents  a 
0.1%  deviation  from  the  OS  curve,  and  the  2”*^  order  correction  deviates  from  the 
order  curve  by  about  0.005%.  It  will  become  evident  in  the  linear  simulations  that  the 
.subtle  differences  between  the  initial  conditions  generated  by  each  of  these  methods 
are  important.  It  should  also  be  noted  that  the  OS  equation  predicts  a  stable  mode 
while  the  corrected  value  represents  an  unstable  one,  for  this  pair  of  a  and  \/^. 
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(a) 
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Deviotion 


Figure  4.2  Initiol  Profiles  for  v.  (o)  Generated  by  OS  analysis. 

(b)  Scaled  deviation  between  1st  order  correction  and 
OS  (dashed),  and  between  2nd  order  correction  and  1st 
order  correction  (dotted). 
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4.2  Linear  Simulations 

The  linear  simulations  without  period  fitting  are  solutions  to  equations  (3.10)  at 
each  time  step.  As  results,  we  present  the  growth  rate,  lm(a;),  as  a  function  of  time, 
where  u  is  computed  at  each  time  step  by  equation  (2.34).  We  note  that  a  positive 
value  of  Im(u>)  indicates  instability.  We  also  point  out  that  if  u  is  fixed  at  its  initial 
profile  during  the  simulation,  and  not  allowed  to  grow,  this  is  equivalent  to  the  parallel 
flow  assumption  and  the  value  of  oj  should  be  constant  in  time  and  equal  to  the  OS 
value.  This  provides  an  initial  check  of  the  simulation  code,  and  is  what  we  observed. 

For  the  simulations  without  period  fitting,  the  dispersive  error  in  the  Crank- 
Nicolson  time  differencing  scheme  makes  the  results  very  sensitive  to  the  number  of 
time  steps  used  per  period  of  oscillation.  We  begin  by  presenting  results  for  a  simu¬ 
lation  which  uses  200  time  steps  per  initial  period  of  oscillation;  t.e..  At  is  chosen  so 
that  based  upon  the  value  of  w  at  t  =  <o,  one  period  of  oscillation  is  divided  into  200 
grid  points: 

“  26oReal(wo)  ' 

For  this  case,  we  give  numerical  solutions  using  each  of  the  three  initial  conditions 
detailed  in  §2.3.  Figure  (4.3a)  gives  the  results  for  the  initial  condition  generated 
by  the  OS  analysis  (solid  line),  emd  the  order  correction.  It  is  evident  that  the 
nonphysical  oscillations  in  the  OS  case  are  a  result  of  the  poor  starting  condition.  The 
1**  order  correction,  which  is  sensitive  to  the  evolution  of  the  mean  flow,  eliminates 
most  of  the  oscillations.  The  simulation  using  the  initial  condition  generated  by  the 
2"**  order  correction  is  free  of  oscillations  on  the  graphiced  sczile.  This  result  is  given 
in  figure  (4.3b). 

Figure  (4.4)  gives  the  solutions  corresponding  to  100  (solid)  and  50  (dashed)  time 
steps  per  initial  period  of  oscillation,  respectively.  The  dotted  curve  on  this  figure  is 
the  200  time  step  per  period  result  given  in  figure  (4.3b).  We  see  that  using  100  time 
steps  per  period  is  still  ample,  but  that  the  quality  of  the  results  when  50  are  used  is 
msu'kedly  diminished  due  to  dispersive  error  in  the  time- differencing  scheme. 
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In  figure  (4.5),  we  give  the  results  for  the  linear  simulations  with  i>eriofl  fitting, 
using  a  time  step  equivalent  to  10  (solid)  and  2  (dashed)  grid  points  per  period, 
respectively.  These  are  .solutions  to  equations  (3.11)  at  each  time  step.  Again,  the 
dotted  curve  ou  this  figure  is  the  200  time  step  per  period  result  of  figure  (4.3b).  We 
note  that  we  get  results  identical  to  this  control  case,  using  as  few  as  2  time  .steps  per 
period. 

We  present  two  final  figures  a,s.sociated  with  the  linear  simulations.  Figure  (4.6) 
shows  Real(u))  versus  time  for  the  simulation  with  period  fitting  using  2  time  steps 
per  period,  and  demonstrates  the  effect  of  the  mean  flow  evolution  on  the  period 
of  oscillation.  In  figure  (4.7)  we  give  a  direct  comparison  between  the  growth  rate 
predicted  by  OS  theory  at  each  time  (da.shed),  versus  the  computed  growth  rate  from 
the  linear  simulation  (solid),  again  using  the  2  time  steps  per  period  case.  We  see  that 
the  evolution  has  a  significant  destabilizing  effect. 
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Figure  4.7  Comporison  of  Orr-Sommerfeld  growth  rote  computed 
locally  at  each  time  step  (dashed)  versus  growth 
rate  computed  from  linear  simulation  which  allows 
growth  of  the  boundary  layer  (solid). 


4.3  Nonlinear  Simulations 

For  all  the  nonlinear  simulations,  we  use  the  initiail  condition  generated  by  the  2'“^ 
order  correction  for  the  j  =  1  mode,  and  its  conjugate  for  the  j  =  —I  mode.  However, 
we  adjust  the  amplitude  to  test  for  the  presence  of  nonlinear  effects.  All  higher  order 
modes  are  initially  taken  to  have  0  values,  as  is  the  mean  flow  correction  uo- 

For  the  first  fovu:  sets  of  simulation  results,  we  used  17  modes  {j  =  —8,— 7,..., 8) 
and  we  fixed  the  peak  amplitude  of  the  initial  profile  for  the  fundamental  mode,  given 
in  figure  (4.2),  to  be  0.000001  (figure  4.8),  0.002  (figure  4.9),  0.005  (figure  4.10)  and 
0.008  (figure  4.11).  We  employed  the  period  fitting  algorithm  and  used  40  time  steps 
per  one  period  of  the  fundamental  mode.  At  each  time  step,  the  set  of  nonlinear 
equations  given  by  (3.18)  and  (3.19),  with  boundary  conditions  (3.16)  and  (3.17)  are 
solved.  Each  of  these  figures  contains  two  graphs:  the  (a)  part  gives  the  growth  rate 
of  the  fundamental  mode,  and  in  the  (b)  part,  we  present  the  energy  curves  for  modes 
j  =  1  to  _)  =  8  where  we  define  energy  as  ln((y(^|),  and  Ej  is  defined  for  the  j*’*  mode 
by  equation  (2.34).  The  strong  oscillations  exhibited  in  the  higher  order  modes  in 
figures  (4.10)  and  (4.11)  result  from  strongly  nonlinear  behaviour  in  the  initial  stages 
of  the  simulation,  since  all  higher  order  modes  eure  initially  set  to  0.  This  effect  is 
amplified  by  aliasing  errors. 

On  each  of  these  figures,  the  dotted  cxirve  represents  the  results  from  the  linear 
simulation  which  uses  the  same  initial  conditions.  For  the  0.000001  case,  the  linear 
results  are  obtained  over  the  entire  time  range  (dotted  curve  sits  on  the  solution  for 
the  fundamental  mode).  However,  even  when  the  eunplitlide  is  as  small  as  0.002,  the 
nonlinear  effects  can  be  seen  within  just  a  few  periods. 

To  test  the  effects  of  even  higher  order  terms,  we  repeat  the  0.008  case  above,  using 
33  modes  {j  =  — 16,  — 15, ...,  16).  Again  we  use  40  time  steps  per  one  period  of  the 
fundamental  mode.  These  results  are  presented  in  figure  (4.12).  In  the  (b)  peirt  of  this 
figure,  we  again  show  only  the  modes  j  =  1  through  j  =  8.  A  reduction  in  oscillations 
occurs  due  to  the  elimination  of  aliasing  errors;  however,  the  startup  effects  can  still 
be  seen. 
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Finally,  to  verify  that  the  period-fitting  algorithm  and  nonlinear  algebraic  solver  are 
working  successfully  for  these  simulations,  we  include  a  simulation  where  explicit  2'“^ 
order  Adams-Bashforth  temporal  discretization  is  used  to  treat  the  nonlinear  terms, 
and  period  fitting  is  not  used.  At  each  time  step,  we  solve  the  linear  set  of  equations 
(3.12)  through  (3.17).  In  this  case,  400  time  steps  per  one  period  of  the  fundemiental 
period  are  used,  and  we  include  16  modes.  Figure  (4.13)  compares  the  growth  rate  of 
this  simulation  (20  periods  in  duration)  with  the  first  20  periods  of  figure  (4.11).  It  is 
evident  that  these  algorithms  are  performing  effectively. 
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Figure  4.8  Nonlinear  Simulation  with  period  fitting.  Durotion:  50 
periods.  17  modes.  Initial  amplitude  of  fundamental 
mode:  0.000001.  (a)  Growth  rote  of  the  fundomental 

mode,  (b)  Energy  curves.  Solid  curves  in  (b)  correspond 
to  the  first  and  second  harmonics:  dashed  curve  is  the 
mean  flow  correction.  Dotted  curve  in  both  figures  is 
the  linear  result. 
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Figure  4.9  Nonlinear  Simulation  with  period  fitting.  Durotion:  50 
periods.  1 7  modes.  Initial  amplitude  of  fundamental 
mode;  0.002.  (a)  Growth  rote  of  the  fundamental  mode, 

(b)  Energy  curves.  Dashed  curve  in  (b)  is  the  mean 
flow  correction.  Dotted  curve  in  both  figures  is  the 
linear  result. 
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Figure  4.10  Nonlinear  Simulation  with  period  fitting.  Duration:  50 
periods.  17  modes.  Initiol  amplitude  of  fundamental 
mode;  0.005.  (a)  Growth  rote  of  the  fundamentol  mode, 

(b)  Energy  curves.  Ooshed  curve  in  (b)  is  the  mean 
flow  correction.  Dotted  curve  in  both  figures  is  the 
linear  result 
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Figure  4.11  Nonlineor  Simulation  with  period  fitting.  Duration:  50 
periods.  17  modes.  Initiol  omplitude  of  fundamental 
mode:  0.008.  (a)  Growth  rate  of  the  fundomentol  mode, 

(b)  Energy  curves.  Doshed  curve  in  (b)  is  the  mean 
flow  correction.  Dotted  curve  in  both  figures  is  the 
linear  result. 
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Figure  4.12  Nonlinear  Simulation  with  period  fitting.  Durotion;  50 
periods.  33  modes.  Initial  amplitude  of  fundamental 
mode:  0.008.  (a)  Growth  rote  of  the  fundomental  mode, 

(b)  Energy  curves.  Dashed  curve  in  (b)  is  the  mean 
flow  correction. 
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Figure  4.1 3  Nonlineor  Simulation  without  period  fitting.  Explicit 
treotment  of  nonlineor  terms.  Durotion;  20  periods. 

17  modes.  Initial  omplitude  of  fundamental  mode:  0.008. 
Growth  rate  of  the  fundomentol  mode  is  shown.  Dotted 
curve  is  simulation  associated  with  figure  (4.1  1)  - 
obscured. 
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§5  Conclusions 

In  this  article  we  have  endeavoured  to  describe  how  Tollmien-Schlichting  waves 
may  develop  in  a  temporally  evolving  flow.  We  have  shown,  as  is  to  be  expected,  that 
above  a  certain  critical  amplitude  their  presence  will  lead  the  situation  to  undergo 
some  kind  of  transition.  We  have  not  tried  to  describe  the  flow  within  this  regime,  we 
have  pmely  tried  to  demonstrate  a  possible  mechanism  which  may  lead  to  this  state 
of  affairs. 

The  main  theme  of  our  article  is  to  develop  reliable  numerical  techniques  which 
allow  us  to  solve  the  nonlinear  partial  differential  equations  that  can  arise  in  this  kind 
of  problem.  The  methods  that  have  been  used  have  been  drawn  from  the  evolution  of 
so-called  parabolising  numerical  techniques,  both  for  instability  work  (see  Bertolotti 
(1991))  and  for  full  Navier-Stokes  simulations  (El-Hady  (1989)).  In  our  work  we  have 
the  distinct  advantage  that  our  underlying  physical  problem  evolves  with  time  rather 
than  a  downstream  coordinate,  which  renders  the  governing  equations  parabolic  rather 
than  elhptic  as  in  the  spatially  evolving  boimdeiry  layer  case.  This  means  that  it  is 
not  necessary  to  change  the  characteristics  of  the  equations  to  effect  a  cheap  numerical 
solution,  rather  the  methods  can  be  exploited  to  provide  more  accurate  solutions  whilst 
using  far  less  resolution,  for  instance  by  period  fitting.  In  the  PSE  work,  the  spatially 
evolving  wavenumber  is  taken  to  vary  linearly  with  the  downstream  coordinate,  so  that 
the  governing  equations  are  no  longer  elliptic.  If  the  next  order  terms  were  retained, 
the  equations  would  revert  to  their  previous  state,  whereas  in  our  problem  we  are  free 
to  retain  as  high  order  variations  as  we  choose.  The  assumption  we  choose  to  make  is 
that  we  can  fit  a  periodicity  and  the  remainder  of  the  evolution  is  on  a  slow  scale. 

In  deriving  initial  conditions  for  our  cedculations  we  start  with  the  Orr-Sommerfeld 
results,  obtained  by  Otto  (1994).  As  mentioned  earlier  we  are  aware  that  these  results 
are  fatally  flawed,  since  they  do  not  allow  for  the  evolution  of  the  undisturbed  flow, 
but  nevertheless  they  provide  an  initial  guess  for  the  modes’  form  and  frequency.  We 
then  consider  an  eigenproblem  in  which  we  not  only  solve  for  the  mode  and  complex 
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frequency,  but  also  the  temporal  derivative  of  both  quaiitities.  In  the  example  con¬ 
sidered,  this  produces  a  fairly  large  change  (rendering  a  mode  that  Orr-Sommerfeld 
predicts  a-s  stable  to  be  unstable).  This  can  be  extended  so  that  we  ensure  the  second 
derivatives  of  the  quantities  are  correct,  producing  minor  modifications  to  the  modes 
chairacteristics  and  eliminating  just  about  all  the  visible  transients  from  the  growth 
rate  curves  (Figure  4.3).  Unfortunately  one  problem  that  we  have  not  resolved  is  the 
transients  due  to  the  start  of  the  nonlinear  modes;  at  the  present  we  start  our  modes 
with  zero  amplitude  (for  the  non-fundamental).  We  believe  that  this  is  what  may  be 
leading  to  the  oscillations  in  the  energy  curves  in  figures  (4.10b),  (4.11b)  and  (4.12b). 

Our  simulations  for  the  linear  problem  seem  to  produce  the  expected  results;  that 
is,  the  flow  can  support  Tollmien-Schlichting  waves  for  a  finite  time  interval.  In  the 
simulation  associated  with  figure  (4.7),  it  appears  that  the  evolution  of  the  mean  flow 
has  a  significant  destabilising  effect,  but  it  should  be  noted  that  the  OS  cuialysis,  to 
which  we  are  making  comparisons,  is  suspect  when  applied  to  this  evolutionary  flow 
situation:  in  Section  4.1,  we  point  out  that  the  OS  analysis  predicts  a  stable  mode 
where  the  corrected  cases  suggest  instability  (we  appreciate  that  this  mode  is  close  to 
the  neutral  boundary).  It  should  also  be  noted  that  even  changing  the  quantity  which 
represents  the  energy  of  the  mode  can  significantly  alter  the  results,  as  seen  in  Smith 
(1979a).  Perhaps  the  main  message  to  be  taken  from  this  p£irt  of  the  study  is  that  one 
should  consider  the  corrected  Orr-Sommerfeld  equation  rather  than  the  parallel  one 
as  the  extra  work  is  more  than  rewarded  by  the  enheinced  accuracy. 

The  nonlinear  simulations  are  as  would  be  expected,  in  that  above  a  certain  initizil 
amplitude  of  the  fundamental,  after  a  period  of  growth,  the  modes  will  induce  a 
transition.  The  development  of  this  set  of  methods  has  enabled  us  to  study  these 
modes,  and  we  are  currently  attempting  to  take  this  calculation  further,  perhaps  on 
into  the  turbialent  regime.  It  is  likely  that  our  modal  calculations  will  have  more 
success  than  the  spatial  calculations  since  our  equations  are  actually  truly  parabolic. 

We  hope  to  extend  this  work  into  the  compressible  domain,  so  that  a  study  of  the 
unstable  characteristics  in  the  case  of  subsonic,  transonic  and  supersonic  plate  speeds 
may  be  performed,  and  the  noise  radiation  at  high  speeds  predicted. 
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